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Abstract 


An analytic model used to predict energy deposition fluctuations in a micro- 
volume by ions through direct events is improved to include indirect delta ray 
events. The new model can now account for the increase in flux at low lineal 
energy’ when the ions are of very’ high energy. Good agreement is obtained 
between the calculated results and available data for laboratory ion beams. 
Comparison of GCR (galactic cosmic ray) flux between Shuttle TEPC. (tissue 
equivalent proportional counter) flight data and current calculations draws a dif- 
ferent assessment of developmental work required for the GCR transport code 
(HZETRN) than previously concluded. 


Introduction 

As the space program enters an era of extended 
manned space operations, the protection of astronauts 
from galactic cosmic ray (GCR) exposure becomes an 
important issue (refs. 1 and 2). The interaction of high- 
energy heavy ions originating in deep space with 
nuclei of shielding structures and body tissues results 
in energy degradation and nuclear fragmentation. 
These nuclear fragmentations produce secondary- and 
subsequent-generation reaction products that alter the 
elemental and isotopic composition of the transported 
radiation fields. The altered radiation level and energy 
spectrum of each ion species can be best estimated 
by using a computationally efficient GCR transport 
code (HZETRN) (refs. 3, 4, and 5) developed at the 
Langley Research Center. Only with the knowledge of 
radiation fields at specified organ locations of astro- 
nauts, can one start to estimate the risk of their expo- 
sure to space radiations. 

Although the HZETRN code was developed by 
using state-of-the-art nuclear models (refs. 6 and 7) 
that had been reasonably tested with laboratory data 
(ref. 8), an integral validation of the code including 
environmental models, atomic and nuclear physics, 
transport algorithm, and vehicle geometry models was 
needed (ref. 9). In previous studies (refs. 9 and 10), we 
compared our results with measurements made with 
a tissue equivalent proportional counter (TEPC) 
onboard the Space Shuttle for both low (28°) and high 
(57°) inclination orbits. The agreement of predicted 
and measured lineal energy spectra for GCR was 
within 70 percent for the region above 2 keV/pm but 
within a factor of 2.3 underpredicted for the region 
below this value. The predicted lineal energy spectra 
were calculated by using HZETRN and postprocess- 
ing the results with an analytic code that describes 


fluctuations in energy deposition from direct ioniza- 
tion events as the ions penetrate through the TEPC 
detector volume. Although the cause for the underpre- 
diction at low lineal energy is conjectured to be a 
neglect of pions in the HZETRN, indirect delta ray 
events should be considered in the TEPC response 
model before any definite conclusion can be made. 

Monte Carlo simulations have been the traditional 
method to model energy deposition by ions in micro- 
volumes. Although such simulations have proven to 
be valuable, they involve sophisticated computer 
codes and time-consuming scoring techniques, require 
large quantities of input information, and are often 
tedious to interpret the results. For complex radiation 
fields such as GCR, it can be a very time-consuming 
process or virtually impossible to use a Monte Carlo 
method to cover all the species and energy range. The 
analytic method of Xapsos et al. (ref. 1 1 ) used in the 
previous studies (refs. 9 and 10) is suitable for such 
complex radiation fields and will be extended here to 
include indirect delta ray events. In this paper, we 
outline the method that describes energy deposition 
fluctuation for both direct and indirect events, com- 
pare our calculated result with available data for single 
ion and energy, and improve our computation in com- 
parison with the Shuttle TEPC data. 

Analytic Model 

The ability of the earlier analytic model by Xapsos 
et al. (ref. 1 1 ) in predicting TEPC response in complex 
radiation fields results from the use and the observa- 
tion that the energy-loss straggling approximates a 
lognonnal distribution with the parameters given in 
terms of relative variance of the random variables 
involved in the energy deposition process. The derived 
(ref. 11) formulations for relative variance depend 



only on easily obtainable macroscopic quantities such 
as LET (linear energy transfer) and range data. 
Herein, we consider only the single event of ions and 
secondary (direct or indirect) delta rays as appropriate 
for the space environment and assume a wall-less 
micro volume. The method outlined begins with the 
calculation of the energy (or ionization) distribution 
for ion (direct) events, followed by an analogous cal- 
culation for electron (indirect) events and, finally, the 
combination of the two distributions. 

Distribution for Ion Events 

When assuming that the ion loses only a small 
fraction of its energy as it traverses through the target, 
the average energy deposited is given by (ref. 1 1 ) 


gling process. A number of treatments of energy-loss 
straggling distributions are well known, such as those 
of Landau (ref. 12) and Vavlilov (ref. 13). However, 
for a given incident ion, each of these distributions is 
applicable only over a limited target size. The use of 
lognormal distribution by Xapsos et al. (ref. 11) was 
first motivated by the work of Condon and Breit 
(ref. 14), Burke (ref. 15), Wilson, Metting, and 
Paretzke (ref. 16), and an abstract by Lepson (ref. 17) 
stating that classical straggling models for high atomic 
number media are not as accurate as results obtained 
using a two-parameter lognormal distribution. Further- 
more, by virtue of central limit theorem (ref. 18), the 
lognonnal distribution results because with each colli- 
sion the ion loses some random fraction of its energy 
that is proportional to its energy before the collision. 


p. = J7 , T . c. (]\ 

ion ion ion ion v ' 

where Z. (on is the LET of the traversing ion; ,v lo|) , the 
path length through the target; and F 10n , the fraction 
of energy initially deposited which remains within the 
site, that is, the fraction not carried out of the site as 
kinetic energy of secondary electrons. The average 
number of electron-hole pairs produced in the target as 
the ion traverses a path length of s ion is then given by 


As described previously, the standard deviation 
and the mean of the lognormal distribution are related 
to x ion and F ion by the relations (ref. 1 1) 

<4, = to 0 + ViJ (4) 

4„„ - (5) 



where lT 10n is the average energy required to produce 
an electron-hole pair. 

The relative variance of energy deposited for the 
ionization event is given by 

F ion = F str,ion + K F,ion < 3 ) 

where F str j on is the energy-loss straggling contribu- 
tion and V /. jo,, is the contribution of Fano fluctua- 
tions. The latter contribution is included if ionization 
is the process of concern as in the case of TEPC. It is 
omitted if energy deposition is the process of concern. 
These two quantities can be detennined as described 
in reference 10. 

One of the key features of the energy deposition 
distribution is how to describe the energy-loss strag- 


where all the variables are a function of ion path 
length s ion . The normalized probability of producing 
x io]1 electron-hole pairs within the target site when the 
ion takes a path length .v ]()n through the site is 


/W-Wfron) 


'^^ x ion°ion 

x ex p{ [ ln (X i°n) “ ^i°n ]2 } (6) 


This equation must be integrated over all possible ion 
path lengths as 

/ion( X ion) = J ^ion C (W *km ^ 

where c(x ion ) is the normalized path length distri- 
bution of the target microvolume. The quantity 
/jo n (X|on) is the normalized probability density that a 
single ion produces x ion electron-hole pairs within the 
restricted target volume upon crossing the site. 
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Distribution for Electron Events 

If the ion misses the target site, there is still a 
probability that energy can be deposited in it by an 
electron (indirect) event. Compared with ion events, 
electron events are more complicated to deal with 
because they have a distribution of energies through- 
out the irradiated medium. The average energy and 
average number of electron-hole pairs deposited in the 
target by an electron traveling a distance s e are given 
by 


£ e = Ve 


W n 


( 8 ) 

( 9 ) 


where L e is the average, slowed electron LET, which 
is obtained by assuming a 1/Z, slo wing-down spec- 
trum and by assuming a ME 1 spectrum with E being 
the electron energy initially produced by the incident 
ion (ref. 19), and W e is the average energy required to 
produce an electron-hole pair. 

The path-length-dependent relative variance of 
ionization for electron events is given by 


y e = 


Vr + V „ +V r 

L,e str,e F,e 


( 10 ) 


from micrometers to nanometers. The smaller the site 
size, the greater is the relative importance of energy- 
loss straggling. Therefore the conclusion can be made 
that the electron ionization distributions are controlled 
by energy -loss straggling and can therefore be approx- 
imated by a lognormal distribution over the complete 
range of interest. We can thus proceed in an analogous 
manner to that in the section “Distribution for Ion 
Events” and obtain equations similar to equations (4) 
through (7) for electron events. 

Combined Distribution 

The remaining problem is how to combine the 
normalized probability density for ion events / ion (x) 
with the normalized probability density for electron 
events / (x) with x being the number of electron- 
hole pairs produced without distinguishing the events. 
The combined ionization distribution / (x) is given by 

/(*) = P fnnbkJ + tt-P) fe(x e ) HD 

where P is the fraction of ion events. We know from 
an earlier work that the fraction of electron events is 
given by (ref. 1 9) 


1 -P = 


O-^ion)* 


( 12 ) 


where V L e is the relative variance of the LET 
distribution; V sUe , relative variance of energy-loss 
straggling; and V F e , relative variance of Fano fluctu- 
ations. Again, these quantities can be calculated as 
described in reference 19. 

Given the complication that the electrons are not 
monoenergetic, it might at first seem as though the 
lognormal distribution would not be directly applica- 
ble. However, examination of figure 7 in reference 15 
indicates otherwise. It clearly shows that the lognor- 
mal distribution closely describes the ionization pro- 
duced by a Co 60 source in a simulated 0.92-pm- 
diameter sphere of tissue-equivalent material. Because 
we now know that the straggling process is described 
by the lognormal distribution, this implies that the 
path length and LET distributions of the Compton 
electrons and photoelectrons are relatively unimpor- 
tant in determining the ionization distribution. In the 
current work, we are concerned with site sizes ranging 


where x is the average number of ionizations of 
combined distribution. Further x can be expressed in 
terms of known quantities as follows: 

1 _ ^ion y ion (13) 

x x ion x e 

where x ion and x e are calculated for the average path 
lengths ,v jon and s e . Equation (1 1) is now ready to be 
tested for single ion and energy before being used in 
predicting Shuttle TEPC results. 

Results 

Monoenergetic Beam 

The earlier model which treats direct events only 
has been tested against experiments and Monte Carlo 
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results for various cases of monoenergetic beam as 
described in reference 10. In general, the model was 
successful in predicting the results except for ions with 
very high energy. Some discrepancy was seen 
between the model and the experiment for a wall-less, 
2-pm-diameter spherical TEPC which measures dose 
distribution in a water column irradiated by 3.9-GeV 
nitrogen ion beam (ref. 20). The underprediction in the 
low y (lineal energy) region is mainly a result of 
neglecting the indirect events since their contribution 
tends to become significant as the beam energy per 
amu increases for a given sensitive volume. For com- 
parison purpose, the earlier result for the case is repro- 
duced here from reference 10. (See fig. 1(a).) The 
current model with both direct and indirect events is 
seen to improve the results greatly. (See fig. 1(b).) 

For a moderately low-energy ion beam, the indi- 
rect event also becomes important when the sensitive 
site decreases to a submicron size. Shown in fig- 
ure 2(a) is the probability density distribution of ion 
pairs produced by 0.3-MeV/amu alpha particles inci- 
dent on a 20-nm-diameter spherical site as calculated 
by Olko and Booz (ref. 2 1 ) by using Monte Carlo sim- 
ulation and also by the current analytic model. The 
model compares reasonably well with the Monte Carlo 
results but seems to slightly underpredict the indirect 
events, and the peak associated with direct events is at 
slightly higher ionization value than the Monte Carlo 
results. In this calculation, the fraction of energy 
deposited in the volume due to ion events F jon is 
about 94 percent. It is notable that the expression for 
/•’ ion as given in reference 22 slightly overestimates in 
the range just below unity when compared with the 
results of Olko and Booz as shown in figures 3 and 4 
of reference 22. It is also shown in reference 21 that 
this fraction varies somewhat among the results 
obtained by Olko and Booz (ref. 21), Kellerer 
(ref. 23), and Guenter and Schultz (ref. 24). Recalcu- 
lation of the case for 0.3-MeV/amu alpha particles is 
thus made by arbitrarily decreasing the fraction to 
92 percent and a much better agreement is reached as 
shown in figure 2(b). Considering the simplicity and 
crudeness of the approach used in evaluating the 
fraction (ref. 22), the slight disagreement shown in fig- 
ure 2(a) is quite reasonable. 


Flight Experiment 

The recalculated lineal energy spectra for GCR 
flux are shown in figures 3(a) and 3(b) for STS-56 
(57° inclination orbit) and STS-51 (28° inclination 
orbit), respectively. These results are obtained by 
using the same version of HZETRN as in previous 
studies (refs. 9 and 1 0) so that the impact of indirect 
events added to the response function can be identified 
when compared with the previous results which are 
reproduced here in figure 4. The improved TEPC 
model tends to raise the predicted flux in the lower y 
region and to lower in the higher y region as was 
expected but still with good agreement for the STS-51 
flight. For the STS-56 flight, there is an overall under- 
estimate across the entire spectrum by the improved 
model. This underestimate suggests there might be a 
systematic problem such as the crudeness in geomag- 
netic transmission function (a simple vertical cutoff 
model) used in our LEO (lower Earth orbit) GCR 
environment of the transport code that could affect 
high inclination orbit more than low inclination orbit. 
The rather large underestimate below 2 keV/pm in 
previous results is now reduced to a level consistent 
throughout the entire spectrum. An improvement to 
the current geomagnetic transmission function is 
needed. 

Concluding Remarks 

The improved analytic model for predicting 
energy deposition fluctuations has been shown to 
agree well with laboratory data and Monte Carlo cal- 
culations for monoenergetic beam even when there is 
a significant contribution from indirect delta ray 
events. These events are important when the ions are 
of very high energy relative to the sensitive site size. 
Comparison of current prediction for GCR (galactic 
cosmic ray) flux in LEO (low Earth orbit) with mea- 
surements obtained by Shuttle TEPC indicates there 
may be an underestimate by the transport code for 
high inclination orbit; this underestimate may result 
from the crudeness of the geomagnetic transmission 
function used in the environmental model and possibly 
other physical processes yet to be included. Still good 
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agreement exists for the low inclination Shuttle flight, 
and the most noticeable difference between the data 
and calculations below 2 keV/|im found in previous 
studies has decreased with the current model. 
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(a) Direct ionization only. (b) Improved analytic model. 

Figure 1. Dose distributions measured by 2-pm-diameter spherical, wall-less TEPC 2-cm depth in water column irradiated by 
3.9-GeV nitrogen ion beam. 



(a)F ion =0.94. (b) -f’ion = 0.92. 

Figure 2. Probability density distribution of ionizations produced by 0.3-MeV/amu alpha particles on a 20-nm-diameter water 
sphere. 
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(a) STS-56. (b) STS-51. 

Figure 3. Calculated GCR spectrum and Shuttle TEPC flight data with TEPC response model used in postprocessing 
HZETRN results including both direct and indirect events. 



10 1 10° 10 1 10 2 10 3 KT 1 10° 10 1 10 2 10 3 

Lineal energy, y, keV/pm Lineal energy, y, keV/pm 

(a) STS-56. (b) STS-51. 


Figure 4. Calculated GCR spectrum and Shuttle TEPC flight data with TEPC response model used in postprocessing 
HZETRN results considering only direct events. 
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